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Abstract 

The  precise  determination  of  a  satellite's  position  is  of  prime  importance  in 
satellite  geodesy  and  in  the  use  of  satellites  in  studying  geodynamics.  It  is  not 
possible,  however,  to  obtain  a  precise  estimate  of  the  satellite's  position  if  suf¬ 
ficiently  large  errors  exist  in  the  gravity  model.  In  this  investigation,  the 
application  of  a  dynamic  model  compensation  method  to  the  problem  of  estimating  the 
motion  of  a  near-earth  satellite  in  the  presence  of  gravity  model  errors  is  described. 
Based  on  numerical  results  obtained  in  a  computer  simulation,  it  is  concluded  that 
the  algorithm  will  yield  an  accurate  estimate  of  the  state  in  the  presence  of 
gravity  model  error:,.  The  effects  of  the  observation  type  and  accuracy,  station 
location  and  observation  batch  size  on  the  accuracy  of  the  orbit  determination 
procedure  are  considered.  Finally,  it  is  shown  that  in  addition  to  an  accurate 
estimate  of  the  position  and  velocity»the  method  can  be  used  to  obtain  an  accurate 
estimate  of  the  unmodeled  acceleration  components. 
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The  precise  determination  of  a  satellite's  position  is  of  prime  importance  in 
satellite  geodesy  and  in  the  use  of  satellites  in  studying  geodynamics.  It  is  not 
possible,  however,  to  obtain  a  precise  estimate  of  the  satellite's  position  if 
sufficiently  large  errors  exist  in  the  gravity  model.  In  this  investigation,  the 
application  of  a  dynamic  model  compensation  methal  to  the  problem  of  estimating  the 
motion  of  a  near-earth  satellite  in  the  presence  of  gravity  model  errors  is  described 
MOM  Based  on  numerical  results  obtained  in  a  computer  simulation,  it  is  concluded 
that  the  algorithm  will  yield  an  accurate  estimate  of  the  state  in  the  presence  of 
gravity  model  errors.  The  effects  of  the  observation  type  and  accuracy,  station 
location  and  observation  batch  size  on  the  accuracy  of  the  orbit  determination 
procedure  are  considered.  Finally,  it  is  shown  that  in  addition  to  an  accurate 
estimate  of  the  position  and  velocity,  the  method  can  be  used  to  obtain  an  accurate 
estimate  of  the  unmodeled  acceleration  components. 
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Introduction 

With  the  recent  improvements  and  refinements  in  the  instruments  used  in  the 
observations  of  near-earth  satellites,  measurements  of  range  and  range-rate  with 
precision  of  less  than  one  meter  and  one  centimeter  per  second,  respectively, 
are  a  reality.  Such  precise  observations  have  potential  for  determination  of  the 
satellite's  state,  namely,  position  and  velocity,  to  the  precision  necessary  for 
geodesy  and  geodyn^mics  studies. 

A  determination  of  the  satellite's  state  commensurate  with  the  observation 
accuracy  using  the  classical  batch  processor  (least  squares)  or  the  sequential 
filtering  method  is  dependent  on  the  accuracy  of  the  tracking  station  locations 
and  the  model  of  the  dynamical  system.  Such  dependence  on  the  model  of  the 
dynamical  system  is  demonstrated  by  Smith,  Kolenkiewicz,  and  Dunn  (1971),  Marsh 
and  Douglas  (1971),  and  Yionoulis,  Heuring,  and  Guier  (1972).  Uncertainties  in 
the  geopotential,  particularly  resonant  coefficients  in  some  cases,  are  an  im¬ 
portant  error  source.  Additional  errors  often  occur  in  the  atmospheric  resist¬ 
ance  modeling.  The  discrepancy  between  the  true  motion  and  the  modeled  motion 
often  limit  the  usefulness  of  the  orbit  determination  results,  both  from  the  point 
of  view  of  using  the  estimates  for  navigation  purposes  as  well  as  using  them  for 
scientific  purposes. 

The  effect  of  gravity  model  errors  on  the  accuracy  of  batch-type  orbit  de¬ 
termination  methods  was  observed  also  in  processing  the  Lunar  Orbiter  range-rate 
observations.  Observation  residuals  were  obtained  which  clearly  showed  a 
periodic  variation  with  an  amplitude  larger  than  the  observation  noise.  The  work 
of  Muller  and  Sjogren  (1971)  and  Sjogren  and  Muller  (1971)  resulted  in  a  nascon 
representation  for  the  gravity  anomalies. 

A  somewhat  different  approach  was  taken  by  Tapley  and  Ingram  (1971)  and 
Ingram  and  Tapley  (1971)  using  a  sequential  estimation  procedure  which  compensated 
fc?  the  model  errors.  In  this  approach,  referred  tc  as  the  Dynamic  Model  Compensa 


(DMC)  method,  the  model  errors  are  assumed  to  consist  of  a  time  correlated 
component  and  a  purely  random  component  and  are  approximated  by  a  first-order 
Gauss-Marko\  process.  Application  of  the  DMC  method  to  the  lunar  orbit  phase 
of  the  Apollo  10  and  11  spacecraft  showed  that  (1)  the  DMC  method  obtained  ob¬ 
servation  residuals  within  the  observation  noise  level  and  (2)  provided  estimates 
of  the  acceleration  due  to  the  model  error  which  could  he  correlated  with  the 
surface  mans  anomalies  described  by  Sjogren  and  Muller  (1071).  Additional  ap¬ 
plications  of  the  DMC  method  have  been  made  by  Tapley  and  Schutz  (1972)  to  the 
problem  of  estimating  unmodeled  forces  due  to  errors  in  the  lunar  potential. 

In  addition,  Connolly  (1972)  considered  application  of  the  method  to  the  orbit 
determination  problem  of  a  near-earth  satellite. 

Other  applications  of  statistical  models  have  been  explored  by  Shaw,  Paul, 
and  Henrikson  (1969)  for  use  in  inertial  navigation  systems.  Shaw  et^  al.  used 
statistical  models  for  the  vertical  deflection  from  gravity  anomaly  models.  More 
recently,  Kasper  (1971)  and  Jordan  (1972)  have  utilized  second-order  Gauss-Markov 
processes  to  represent  the  gravity  anomaly. 

In  this  investigation,  the  question  of  how  accurately  both  the  state  of  the 
satellite  and  the  acceleration  due  to  model  errors  can  be  estimated  for  a  near- 
earth  satellite  is  considered.  A  near-earth  orbit  of  1000  km  inclined  at  ^2°  is 
assumed  for  the  study.  In  the  investigation,  computer  simulation  techniques  were 
used  to  investigate  the  behavior  of  the  DMC  method  as  applied  to  the  orbit  de¬ 
termination  problem  of  a  near-earth  satellite.  In  this  study,  the  model  errors 
were  assumed  to  arise  from  an  incomplete  description  of  the  geopotential.  Two 
sets  of  simulated  observations  were  generated,  each  resulting  from  a  different 
method  of  imulat  in.T  gravity  model  errors.  It  is  shown  that  state  estimates  cm 

be  obtained  which  reduce  the  observation  residuals  to  within  the  apriori  observa- 

• 

tion  accuracy  and  that  the  methods  can  be  used  to  obtain  an  accurate  estimate  of 
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tin.*  accelerations  due  to  gravity  model  error.  The  influence  oT  observation  ac¬ 
curacy,  the  magnitude  of  the  model  error,  and  the  number  of  tracking  station:;  on 
the  accuracy  with  which  the  acceleration  due  to  the  model  error  can  be  estimated 
in  discussed. 

The  tls I  i di.it  i on  rponMluri' 

The  equation:;  which  describe  the  motion  of  an  earth  satellite  can 
be  expressed  by  the  following  system  of  first-order  equations: 


r  =  v  ,  v  =  a_  t  a_  +  m(t) 
c  p 


(1) 


where  r  is  a  three-vector  of  position  components ,  v  is  a  three-vector  of  vel¬ 


ocity  components,  a£  is  the  acceleration  due  to  the  central  body,  a^  is  the  modeled 

acceleration  due  to  other  sources,  gravitational  or  otherwise,  and  m  represents 

the  acceleration  due  to  incorrect  modeling  in  a  and  a  .  The  term  m(t)  will  be 

c  p 

referred  to  as  the  "unmodeled  acceleration". 

The  three-vector  m(t)  represents  the  effects  of  all  accelerations  not  ;c- 
counted  for  in  the  mathematical  model  used  to  describe  the  motion  of  the  satellite. 
In  this  discussion,  m(t)  is  approximated  as  a  first-order  Gauss-Markov  process, 
c(t),  which  satisfies  the  following  vector  differential  equation 

c(t)  =  Be(t)  +  u(t)  (2) 

where  u(t)  is  a  three-vector  of  Gaussian  noise  whose  components  are  assumed  to 
be  described  by  the  statistics: 


E[u(t)]  =  0  ,  E[u(t)u  (t)]  =  q(t)6(t-x) 

where  q(t)  is  a  3  x  3  positive  definite  matrix.  The  coefficient  matrix  B  is 


(3) 


defined  by  the  components  B..  =  where  T.  ,  the  time- correlation  coefficients, 

7  1  i]  T.  it  l 

l 

are  assumed  to  be  unknown  parameters  whose  values  are  to  be  determined  during  the 
estimation  process.  The  state  vector  is  then  defined  to  be 


VT  r  T  T 
X  =  [r  :  v 


T  _T, 
e  :  T  ] 
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Then  the  differential  equations  of  state  become 

v 

a  ♦  a  ♦  e 
c  p 

Be  +•  u 

0 

with  the  initial  conditions  X(t  )  =  X  . 

o  o 

In  the  following  discussion,  it  is  assumed  that  the  observation  noise,  v^, 
satisfies  the  following  conditions: 

E[vi]  =  0  ,  E[v.v*]  =  Ri«.j  ,  E[v.xj]  =  0  ,  (5) 

where  E[  ]  is  the  expected  value  operator.  Further,  the  observations  are  re¬ 
lated  to  the  state  by  the  relation 


y.  =  G(xi,t.)  +  v.  (6 

where  <»(Xj,t)  i:j  some  nonlinear  function  of  the  state. 

The  problem  considered  then  is  posed  as  follows:  Given  the  relation  for 
propagating  the  state,  Eq.  (4);  the  observation  state  relation,  Eq.  (6);  a 
sequence  of  observations  Y^,  i  =  1,  ...,  k;  the  statistics  on  the  state  noise, 
Eq.  (3);  and  the  observation  noise,  Eq.  (5),  determine  the  best  estimate  of  the 
state  in  the  minimum  variance  sense,  X^,  at  the  time  t^.  Under  these  conditions 
the  estimate  can  be  obtained  using  the  following  algorithm  (Jazwinski,  1969): 

*k  =  9(Xk-l’  Vl’  tk) 

Pk  =  *(tk»  ^^k*  tk  i'  +  \  i 

•Sc  =  K  *  V* 

G(V  V1 


pk  =  [I  -  W  pk  ( 

where  Hk  =  [3G/3XJc]*,  \  =  Ep^  |  Y^  ...,  Y^]  and  \  =  *l\  \  Yj. . Y^.  J 
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The  covariance  matrices  P.  and  P,  are  associated  with  the  state  estimates  X,  and 

k  k  k 

A 

X^,  respectively.  The  state  transition  matrix,  <f>( tk ,tkl)  satisfies  the  follow- 
in)’  differential  equation 

♦(t,  tk)  =  A(t )  ♦(t,  tR)  ;  4(t,  tk)  =  1  (8) 

where  A(t)  =  [3F/3X]*.  The  symbol  []*  indicates  tfiat  the  quantity  in  the  bracket 

A 

is  evaluated  on  the  nominal  solution,  X(t)  =  0(Xk  t)»  *  t,  that  is, 

A 

the  solution  obtained  by  integrating  Eqs.  (4)  with  the  condition  x*(tk  ^)  =  X^^ 
and  u*  =  E[u(t)]  =  0.  The  algorithm  given  by  Eqs.  (7)  and  (8)  is  the  extended 
form  of  the  linear  sequential  estimator,  or  the  Kalman-Bucy  filter,  as  iiscussed 
by  Jazwinski  (1969).  The  covariance  matrix  Pk  is  associated  with  the  best  estimate 
of  Xj^  based  on  k-observations  while  Pk  is  the  covariance  matrix  associated  with 
the  best  estimate  of  Xj^  based  on  (k-1)  observations. 

A 

It  should  be  recalled  that  the  estimate  ,  includes  an  estimate  at  the  time, 

tk,  of  the  components  of  the  position,  the  velocity,  the  unmodeled  acceleration 

c(t)  and  the  correlation  coefficients,  T  ,  T  ,  and  T  .  The  algorithm  requires 

x  y  z 

apriori  or  "initial"  estimates  of  each  of  these  quantities  as  well  as  the  apriori 

covariance  matrices,  P  ,  Q.,  and  R. ,  associated  with  the  initial  state  and  with 

o  1  i 

the  state  and  the  observation  noise,  respectively.  The  development  of  the  algo¬ 
rithm  and  the  computational  procedure  required  to  implement  the  algorithm  are 
discussed  in  greater  detail  by  Ingram  (1971). 

The  following  section  describes  the  generation  of  the  observations  used  in 
the  estimation  procedure  and  discusses  the  simulation  of  the  gravity  model  errors. 
Generation  of  the  Simulated  Observations 

The  simulated  observations  are  generated  by  the  integration  of  the  system 


of  equations 


e 


r\  =  v_  ,  v_  =  a  +  m  (9) 

T  T  *  T  m 

where  r^,  is  a  three-vector  of  the. true  position  components,  v^,  is  a  three-vector 

of  the  true  velocity  components,  a^  represents  the  terms  modeled  in  the  estimation 

algorithm,  arid  m  represents  the  terms  not  directly  included  in  tlio  estimation 

differential  equations.  The  term  m  represents  a  model  error  in  the  estimation 

procedure.  To  simulate  model  errors,  two  forms  for  m  were  used.  First,  m  is 

assumed  to  represent  the  acceleration  due  to  a  point  mass  on  the  surface  of  the 

earth  and  close  to  the  satellite's  ground  track.  The  term  a  includes  all  of  the 

m 

zonal  harmonics  up  to  the  sixth  degree  as  well  as  the  tesseral  and  sectorial 

harmonics  up  to  fourth-order  and  the  lunar  gravitational  effect.  In  the  second 

form,  a  includes  J0  and  the  lunar  effect ,  whereas  m  contained  the  geopotential 
m  2 

terms  up  to  sixth  degree  in  the  zonal  harmonics  and  all  other  terms  to  the  fourth 
order. 

Integration  of  Cq.  (9)  yields  a  solution  for  the  motion  of  the  satellite 
used  in  generating  observations.  The  computer  program  ascertains  which  tracking 
stations  can  observe  the  satellite  and  an  observation  (either  range,  azimuth, 
elevation,  range-rate  or  any  combination)  is  computed.  These  observations  are 
then  corrupted  by  adding  the  quantity,  oA,  where  o  is  the  standard  deviation  of 
the  observation  noise  and  A  is  a  Gaussian  distributed  random ‘variable  with  zero 
mean. 

The  numerical  integration  is  performed  using  the  Runge-Kutta  methods 
formulated  by  Fehlberg  (1968).  As  implemented  in  the  simulation  program,  two  dif¬ 
ferent  order  integrators  can  be  used.  For  periods  when  the  satellite  cannot  be 
seen  by  any  of  the  tracking  stations,  the  Runge-Kutta-Fehlberg  7(8)  is  used  with 
a  variable  stepsize  control.  During  periods  in  which  observations  are  being 
generated,  a  fixed  step  is  used  to  comply  with  a  fixed  observation  interval. 
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This  interval  is  normally  small  enough  to  allow  use  of  a  lower  order  integrator, 
namely,  the  Runge-Kutta-Fchlberg  4(5).  If  the  observation  interval  results  in 
a  stepsize  which  is  too  large  and  the  4(5)  is  inadequate  for  maintaining  accuracy, 
the  7(8)  method  can  be  used. 

Data 

The  estimation  program  requires  the  numerical  integration  of  Eq.  (4)  and 
Eq.  (8),  whereas  the  observation  generation  program  requires  the  integration  of 
Eq.  (9).  The  state  vector  of  Eq.  (4)  consists  of  12  elements,  the  state  transi¬ 
tion  matrix  consists  of  144  elements,  and  Eq.  (8)  consists  of  six  elements. 

Eq.  (8),  can,  however,  be  reduced  to  a  lower  order  by  noting  that  a  number  of 
elements  arc  zero. 

The  initial  conditions  for  Eq.  (4)  and  Eq.  (9)  used  in  the  simulation  are 
shown  in  Table  I.  A  rectangular  coordinate  system  is  used  in  which  X  and  Y  are 
in  the  equatorial  plane  and  Z  is  along  the  polar  axis.  The  ground  track  of  the 
satellite  is  shown  in  Fig.  1  where  the  semi-major  axis  of  the  satellite  orbit  is 
7,507,000  meters,  the  eccentricity  is  .025,  and  the  inclination  is  41.2°.  These 
orbital  elements  closely  approximate  those  of  the  Beacon  Explorer  C  (BE-C) 
spacecraft  (Kolenkiewicz,  1972). 

The  assumed  tracking  station  locations  used  in  the  simulations  are  shown  in 
Table  II.  In  general,  the  observations  from  horizon-to-horizon  of  the  tracking  station 
are  used  in  the  estimation  process.  However,  several  cases  were  examined  in  which  a 
25°  elevation  constraint  was  imposed.  These  results  are  discussed  in  a  later  section. 
The  period  of  time  during  which  each  station  can  observe  the  satellite  is  shown 
also  in  Figure  II.  Although  Figure  II  shows  15  stations,  not  all  stations  are  used. 

The  geopotential  coefficients  used  were  given  by  Melbourne,  et_  al_  (1968).  Since 
these  coefficients  were  being  used  to  simulate  a  gravity  model  error,  a  more  ac¬ 
curate  set  such  as  those  determined  by  Gaposchkin  and  Lambeck  (1970)  was  not  required. 
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In  all  cases ,  the  initial  covariance  matrix,  P  ,  had  the  following  forr. 


r  o2 


r  o2 


r  o2 


v  o2 


Vo2 

* 

°T 

T 


wlieiv  (i  =  1000  meter::,  »i  -  1  meter/sec,  a,.  =  .Oj.5  seconds ,  u,  =  .0001  meters/r.oc2. 
r  v  r  f. 

In  general,  the  initial  value  of  £  is  taken  to  be  zero  unless  otherwise  stated. 

Results 

Not  only  were  two  types  of  gravity  model  errors  studied,  but  the  effects  of  ob¬ 
servation  noise,  station  location,  and  number  of  stations  were  investigated.  These 
results  are  discussed  in  the  following.  The  gravity  model  error  due  to  a  point  mass 
will  be  referred  to  as  Case  I  and  the  second  type  discussed  in  preceding  sections  will 
be  referred  to  as  Case  II.  Furthermore,  Case  I  is  only  considered  in  the  50  cm  range 
observations  study. 

1.  50-cm  Range  Observations 

Both  types  of  gravity  model  errors  were  used  in  the  study  with  range  observa¬ 
tions  in  the  noise  level  of  50  cm. 

A.  Case  I 

In  the  generation  of  the  observations,  a  point  mass  of  3.64  x  1C-5  times 
the  mass  of  the  earth  was  placed  at  41°N  latitude  and  108°W  longitude,  ihe 
standard  deviation  in  the  state  noise  was  .005  and  the  initial  value  assumed  ^er 
the  time  correlation  coefficient,  T^,  i  =  1,  2,  3,  was  200  seconds.  Figure  2 
shows  the  result  of  estimating  the  acceleration  due  to  the  model  error,  or 


uiiiihmIc  1  I  .iccclrrat  ion.  In  this  figutv,  tin;  smooth  1  i lit;  represents  the  true 
value  and  the  jagged  line  represents  the  estimate,  figs.  2a,  2b,  and  2c 
represent  the  acceleration  components  in  a  rectangular  coordinate  system. 

Note  that  the  estimate  follows  the  true  acceleration  component  and  that  the 
estimate  has  a  mean  which  is  close  to  the  true  value.  Figs.  2d  and  2e,  show 
the  trace  of  the  state  error  covariance  matrix  and  the  true  error  norm  plotted 
as  a  function  of  time.  In  Fig.  2d,  for  example,  the  smooth  line  represents 
the  trace  of  P  associated  with  the  position  components  and  the  jagged  line 
represents  the  true  error  norm.  It  is  important  tc  note  that  the  state  error 
covariance  matrix  bounds  the  true  error.  In  this  case,  1?  tracking  stations 
were  used,  although  not  all  stations  can  observe  the  satellite  simultaneously 
as  shown  in  Table  II.  Those  not  used  were  Merritt  Island,  Arecibo,  and 
Bermuda.  Since  the  true  error  is  within  one  meter  from  5  to  20  minues ,  the 
(O-C)’s  would  also  be  within  one  meter. 

B.  Case  II 

The  gravity  model  error  was  simulated  using  J_  through  J  and  the  tesseral 

3  b 

and  sectorial  harmonics  through  order  four.  The  same  stations  used  in  Case  I 
were  also  used  in  this  case.  Figure  3a,  3b,  and  3c  show  the  acceleration 
er.tim.Ho  un-.I  the  true  acceleration.  The  standard  deviation  in  the  state  noise 
was  0.00002b.  Again,  the  smooth  line  is  the  true  and  the  jagged  is  the 
estimate.  Although  the  estimate  of  the  acceleration  is  not  as  good  as  in 
Case  I,  the  state  estimate  is  quite  good,  that  is,  within  one  meter.  It 
should  be  noted  that  the  acceleration  error  is  approximately  100  times 
smaller  than  ir  the  preceding  case.  The  acceleration  error  in  Case  II  is  cn 
the  order  of  5-8  milligals,  as  can  be  seen  in  Fig.  3.  The  (O-C)'s  are  shown 
in  Fig.  4.  Note  that  most  of  the  points  are  within  the  50  cm  noise  level  of 


observation. 
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As  a  point  of  comparison  the  sane  case  was  run  using,  the  hatch  processor. 
The  (O-C)'r,  are  shown  in  fig.  5.  Note  the  periodic-type  behavior  with  ampli¬ 
tude  conn idcruh iy  .ihovo  the  noise  level  of  the  observation. 

Figure  fi  shows  the  results  obtained  by  replacing  the  fhinta  Barbara,  Las 
Vega.-.,  and  Albuquerque  stations  by  ones  at  Merritt  Island,  Arecibo,  and  Bermuda 
The  influence  of  additional  east  coast  stations  was  not  significant  except 
toward  the  end  of  the  pass.  It  can  be  seen  in  Fig.  2d  that  the  error  begins 
a  significant  increase  at  approximately  20  minutes.  Table  II  shows  that  at 
about  20  minutes  there  are  four  stations  taking  data  and  beginning  at  24 
minutes,  only  one  station  taking  data.  The  error  growth  appears  to  be  an  ef¬ 
fect  of  the  number  of  stations.  This  effect  will  be  discussed  in  a  later 
section. 

2.  1-em  Rang#*  Observations 

Repeating  the  preceding  Case  II  with  o  =  50  cm  replaced  by  a  =1  cm  yielded 
the  result:;  shown  in  Tig.  7.  The  tracking  stations  not  used  were  Merritt  Island, 
Arecibo,  and  Bermuda.  It  is  apparent  that  Dy  improving  the  observation  accuracy, 
the  estinate  of  the  unmodeled  acceleration  is  improved  also.  The  upward  trend  of 
the  covariance  trace  in  Fig.  7d  is  again  an  effect  of  a  single  station  observing 
the  satellite.  Note,  however,  that  the  true  error  is  no*  increasing.  ”ic\ce  3  il¬ 
lustrates  the  (O-C)’s  for  two  of  the  stations,  Goldstone  and  Goddard.  V.zzs  that 
most  of  the  residuals  are  within  the  1-cn  observation  noise,  for  comparison,  the 
same  case  was  run  in  a  batch*  processor.  The  (O-C)'s  resulting  from  this  case  are 
shown  in  Fig.  0. 

Since  it  is  evident  that  the  accuracy  with  which  the  unmodeled  acceleration 
can  be  estimated  is  dependent  on  the  observation  accuracy.  Case  I  was  not  run  with 
1— cm  observat ions  since  50— cm  observations  gave  good  estimates.  The  1— cm  ebservati 
would  simply  reduce  the  scatter  in  the  estimate. 
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3.  2mn./t>ec  Range-Rate  Observations 

The  Case  II  discussed  in  thq  preceding  section  (2)  was  repeated  using  range- 

rate  observations  of  2  nsn/sec  accuracy.  The  results  are  shown  in  Fig.  10.  The 

accuracy  with  which  the  unmodelcd  acceleration  can  be  estimated  is  comparable  to  the 
J-cnt  range  observ.il  ion::.  The  (O-C)'s  behave  as  expected. 

4.  25°  elevation  Constraint 

A  case  was  run  to  examine  the  effect  of  a  25°  elevation  constraint  on  the  1-cm 
range  observations.  All  observations  yielding  an  elevation  less  than  25°  were  ’dis¬ 
carded.  The  observation  interval  us^  i  was  .02  min  =  1.2  sec.  Again,  Merritt  Island, 
Arecibo,  and  Bermuda  were  not  included.  In  this  case,  as  can  be  seen  in  Table  II, 
the  first  station  which  can  see  the  satellite  is  San  Francisco  at  approximately  3.6 

minutes.  Shortly  thereafter,  Santa  Barbara  acquires  the  satellite.  Figure  11a  shows 

that  the  acceleration  estimate  is  poor  until  approximately  4  minutes  at  which  time  a 
third  station  acquires,  namely,  Ooldr.tonc.  After  thin,  the  estimator  yield;*  good 
estimate:;.  Thu:,  for  the  simulation,  the  horizon-to-horizon  obr.ervat  ions  yield 
results  representative  of  a  more  realistic  case. 

5.  Effect  of  Number  of  Stations  (..-serving 

The  preceding  case  indicated  the  effect  of  the  number  of  tracking  stations  which 
can  observe  the  satellite.  Figure  12  shows  the  results  of  using  only  thre<  stations, 
namely,  San  Francisco,  Ft.  Davis,  and  Goddard.  The  observation  accuracy  is  1-cm  in 
range,  and  the  observations  are  taken  from  horizon-to-horizon.  Note  the  radical  change 
in  the  Y-component  (Fig.  12b)  when  the  third  station  is  included  at  approximately  8 
minutes.  Even  after  the  San  Francisco  station  is  lost  at  16  minutes,  the  estimate  re¬ 
main-:  quite  good  until  the  Ft.  Davis  station  loses  acquisition.  Thus,  this  result 
shows  that  three  stations  can  yield  a  good  result  while  two  stations  nay  yield  com¬ 
parable  or  reasonable  results.  The  question  of  what  a  single  station  can  do  with 
this  method  was  not  thoroughly  investigated. 


Summary 


Based  on  the  numerical  results  in  the  previous  section,  it  can  be  concluded 
that  the  DKC  method  provides  a  significant  improvement  over  the  conventional  batch 
method  in  estimating  the  state  of  a  satellite  in  the  presence  of  gravity  model 
errors.  The  50  cm  range  observations  and  2  min/sec  range-rate  observations  are 
realistic  while  the  1  cm  range  observation  will  probably  be  available  in  the  very 
near  future.  Using  these  observation  accuracies,  it  was  shown  in  the  previous  sec¬ 
tion  that  the  accuracy  with  which  the  unmodeled  aocolerat ion  can  be  estimated  de¬ 
pends  on  the  observation  accuracy  and  the  magnitude  of  the  acceleration  due  to  the 
model  error.  Even  when  good  acceleration  estimates  are  not  obtained,  the  state  estimate 
yields  (O-C)'s  within  the  observation  noise.  The  number  of  tracking  stations  which 
can  observe  the  satellite  influences  the  acceleration  estimate.  The  results  indicate 
that  two  or  three  stations  yield  good  estimates  at  the  1-cm  range  noise  level. 

Further  work  is  in  progress  to  determine  whether  the  second-order  Gauss-Markov 
process  would  yield  a  significant  improvement  over  the  first-order  process  used  in 
this  paper.  Additional  investigation  of  the  number  of  tracking  stations  required  as 
well  as  the  effect  of  tracking  station  location  errors  is  underway.  Consideration  is 
also  being  given  in  the  use  of  the  acceleration  estimates  to  improve  the  mathematical 
model  of  the  geopotential. 
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Table  II 

Simulation  Tracking  Station  Locations  and  Coverage 
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Velocity  Error  Nom  (meters /sec) 


Position  Error  Norm  (meters)  e.  Velocity  Error  Nona  (raeters/sec 
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Figure  4.  (0-C)  For  thel  DMC  Method 
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Figure  12.  Case  II  with  Three  Tracking  Stations 
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